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Abstract 

Quantum random walks have received much interest due to their non-intuitive dynamics, which 
may hold the key to a new generation of quantum algorithms. What remains a major challenge is 
a physical realization that is experimentally viable and not limited to special connectivity criteria. 
We present a scheme for walking on arbitrarily complex graphs, which can be realized using a 
variety of quantum systems such as a BEC trapped inside an optical lattice. This scheme is 
particularly elegant since the walker is not required to physically step between the nodes; only 
flipping coins is sufficient. 
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Random walks have been employed in virtual 
everyday phenomena such as the DNA synapsis 



ly every science related discipline to model 
1|, animals' foraging strategies Q], diffusion 



and nrobUity in material H and exchange rate forecast Q. They have a.so found algorith- 
mic applications, for example, in solving differential equations |5J , quantum monte carlo for 
solving the many body Schrodinger equation [6|, optimization 7|, clustering and classifica- 
tion fractal theory Q or even estimating the relative sizes of Google, MSN and Yahoo 
search engines [l(]]. Whilst the so called classical random walks have been successfully uti- 
lized in such a diverse range of applications, quantum random walks are expected to provide 



us with a new paradigm for solving many practical problems more efficiently 11 



12j. In fact 



quantum walks have already inspired efficient algorithms with applications in connectivity 
and graph theory [3, Q , as well as quantum search and element distinctness Q , due 
to their non-intuitive and markedly different properties including faster mixing and hitting 
times. 

The question we address in this paper is how to physically implement a quantum random 
walk in the laboratory. Over the last few years there have been several proposals for such a 

3] , ion traps 



physical implementation using Nuc 



191 ]. classical and quantum optics 



quantum dots 



23 



ear Magnetic Resonance [17|], cavity Q 



12,1201], optical lattice and microtraps 



ID 
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as well as 



241 ] . None of the existing proposals however consider quantum random 



walks on general graphs, with the majority describing only a one-dimensional implementa- 
tion. This is while from an application point of view most useful algorithms would involve 
traversing graphs with arbitrarily complex structures. 

In this paper, we present a scheme which considerably simplifies the evolution of the 
quantum walk on a general undirected graph. We then describe a systematic procedure 
capable of performing this quantum walk on a variety of existing as well as prospective 
quantum computing platforms. Finally we present an example of one such implementation, 
using a Bose-Einstein condensate (BEC) of 87 Rb atoms trapped inside a 2D optical lattice 
25]. 

First we consider a complete graph with all possible connections between the M nodes 
including self loops (Fig. [Th). Here the walker requires an A/"-sided coin for moving from 
one node to M other nodes. The complete state of the walker is therefore described by 
1^) = Y^ij=i Y^k=i ■Aj,k\j, where Aj,k are complex amplitudes, \j) and \k) represent the 
node and coin states respectively. A quantum coin flip corresponds to a unitary rotation of 
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the coin states at every node j using an M x M matrix dj also known as the coin operator. 
The coin operation is followed by the walker stepping from node j simultaneously to all other 
nodes on the graph using a conditional translation operator T such that T\j, k) — > \j', k'), 
where j and j' label the two nodes at the end of an edge e^y 26] . The quantum walk evolves 
via repeated applications of the coin followed by the translation operator. More explicitly, 
we have \ip n ) = T n C n . . . T 2 C 2 T\ C\ |^o), where |^o) is the initial state of the walker, |-0„) 
is its state after n steps, Ci and are the coin and translation operators at the ith step, 
and C incorporates the individual coin operators c,\ . . . which simultaneously act on all 
the nodes. The operators c can in principle invoke different rotations at each node j, but 
are often uniformly set to be the Hadamard rotation. 

In traversing the edge e^y, we define T\j, k) — > \k,j) (Fig. [Tb). Without undue loss of 
generality, this choice of translation operator has the unique advantage of being independent 
of graph connectivity, and thus enabling a quantum walk to be systematically implemented 
on any arbitrary graph. Upon visualizing the Hilbert space of the walk as an M x M square 
array Ti with entries hjk representing the states \j, k), the application of the translation 
operator T to the state space of the walk simply becomes equivalent to a transposition of 
the array elements. Let us now consider the first few steps in the evolution of a quantum 
walk. Applying C\ to the state space of the walk involves performing M simultaneous unitary 
transformations Cj, each on the coin states of the node corresponding to the j'th row. This 
leads to a natural grouping of the states along the rows of 7i and we employ the relabeled 
operator to highlight that it operates on horizontally grouped states (Fig. [Tb). What is 
particularly convenient now is that instead of transposing 7i due to the action of 71 we can 
simply transpose the application of the next coin operator C 2 . By transposing C 2 we mean 
regrouping the states, this time along the columns of 7i, and performing M simultaneous 
unitary transformations cy each on the states of the jth column. As before we employ the 
relabeled operator C\ to highlight that it operates on vertically grouped states. 

In the above formulation, the effect of the translation operator T is implicit in the re- 
grouping of states and does not appear in the expression governing the evolution of the walk, 
which can now be written as \ip n ) = C^C^_ X . . .C\ C± \ipo), halving the number of required 
operations. It is in this sense that we have qualified this process as a "quantum random 
walk without walking" ; the walker is not required to physically step between the nodes, only 
flipping the coin is sufficient. As we will see, removing the quantum walk's dependence on 
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the translation operator T greatly facilitates its physical implementation. 

We now construct our intended graph Q by simply removing all the unwanted edges 
(dotted lines in Fig. from its complete counterpart {? max - In turn this has the effect of 
removing some of the states from the Hilbert space Ti (dotted circles in Fig. [lb). Removing 
the edge e 3 y for example, corresponds to removing two states and \j',j)- In our 

approach however, instead of removing these unwanted states from 7i, we simply isolate them 
from interaction with other states by appropriately designing the coin operators c\ . . . c^f. 
Taking C H as an example, matrix cf performs a unitary transformation on the jth row of 
H. Hence to isolate the state \ j, k) we obtain a modified coin matrix whose column elements 
C\k ■ ■ ■ CHk and row elements c^x ■ ■ ■ Ckj\f are all set to zero except for Ckk which is 1. Using 
this modified coin matrix guarantees that if initially the walker has no amplitude in state 
\j, k) , this state will remain unpopulated throughout the evolution of the walk. 

It is clear from the preceding discussion that a physical implementation of this walk 
requires two basic properties commonly found in a variety of systems proposed for traditional 
quantum computing: (a) M 2 basis states arranged in a square array formation and (b) 
implementing the operators (cj>), which at once perform an A/"- state unitary rotation on 
all the amplitudes in row j (column j') of the 2D state space. Such a mechanism can indeed 
be efficiently constructed if the system is capable of performing pairwise unitary operations 
on non-neighboring states similar to those demonstrated in 
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l29l . l30f | and discussed 
321 ] and references therein. The key to our implementation is a Cosine Sine (CS) 



decomposition 



33j which effectively takes the single unitary operator (cj,) and replaces 
it with a series of pairwise operators which we know how to implement. One requirement 
of this implementation is that M = 2 for some integer N, which can introduce some 
redundancy in the the Hilbert space of the quantum walk, but only adds a linear overhead. 
Considering the wave function along row j, we represent the operator Cj cLS cLll 
unitary matrix acting on a vector = (ai ■ • -aw) of amplitudes in row j. Performing 
N — 1 recursive CS decompositions on we obtain 



Ulx 



cf = J] Uiidi), where U^d,) = ^ u i)2 j (1) 

and Ui t k represent d% x di square blocks along the lAi diagonal with k — 1, 2 • • ■ M jd{. Block 
dimensions can vary for each Ui with values restricted to di = 2, 4, 8 • ■ ■ M jl. For = 2, 
blocks u it k represent general 2x2 unitary matrices, but for di > 2 they assume the special 



form 



\ 




















(2) 



l . k 



where each quadrant is diagonal with respective entries c r and s r corresponding to cos(0 r ) 
and sin(0 r ) for some angle <f) r and r = 1,2- • -d/2. The action of each matrix Uiidi) on 
the vector can now be directly implemented using pairwise interactions. Upon a closer 
examination of u^k in Eq. [2]we find that each c r s r square block (dotted) performs a pairwise 
unitary transformation Ui^, r on the amplitudes a^k-i)d+r and a^k-i)d+r+d/2, which are non- 
neighboring for d > 2. Hence the rotation Uiid^) can be applied at once by simultaneously 
activating pairwise interactions between all states in the range \j, kd — d + 1) . . . \ j, kd — 
d/2) and their corresponding counterparts \j, kd — d/2 + 1) . . . \j, kd) for all k. Note that 
conveniently, all interacting pairs of states have the same interval d/2 which greatly facilitates 
the design of a physical implementation. 

In the following we describe one such physical implementation using a BEC trapped in 
a 2D optical lattice 3^], where states \j, k) of the walk are encoded using the individual 
trapping sites and the BEC wave function acts as the quantum walker with some initial 
distribution throughout the lattice sites. The system is driven into a Mott insulator phase 
35J thereby suppressing the tunneling between neighboring lattice sites. A series of specially 
tailored control laser operations are then introduced to address, manipulate and interact the 
BEC wave packets in individual sites, in a way that corresponds exactly to the action of the 
operators Cj (cj,) along the lattice rows (columns). Although the control laser wavelength 
and the lattice period Ai at ti C c are comparable in size, problems associated with unwanted 
interactions of the control laser with nei ghb oring sites can be circumvented by adopting 
techniques such as those detailed in 36|, |37|] or more readily by choosing every 2nd, 3rd 
or £th lattice site to represent the walk states. From an application point of view one 
would commonly start with the BEC entirely localized in one site or uniformly loaded into 



every £th site using pattern loading 



38 



or by employing a recently developed imaging and 

3, 3- The desi S n of a11 



manipulation technique based on scanning electron microscopy 
subsequent control operations ensures that the initially empty intermediate sites would, 
in principle, remain unpopulated throughout the walk. In practice however the spatial 
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separation t also acts as a buffer zone to contain any spiling of the BEC out of its confinement 
lattice-site due to unavoidable experimental imperfections. 

To manipulate the trapped BEC wave packet at a given lattice site we propose performing 
arbitrary unitary transformations on the internal states |0) = \F — I, nip = 1) and |1) = 
\F = 2, nip = 2) of the BEC with the aid of a pair of three-photon Stimulated Raman 
Adiabatic Passage (STIRAP) operations 41] . Each STIRAP requires the use of three control 
lasers (with wavelengths ~ Ai at ticc) applied in the counter intuitive order to transfer the 
atomic population in states |0) and 1 1) , to and from an auxiliary state \a) = \F = 2, nip = 0), 
via an intermediate upper state \u) = \F' = l,nip = 1) that does not get populated during 
the transfer (Fig. [2]). The two-photon A STIRAP |1) < — > \u) < — > \a) has already been 
experimentally demonstrated using circularly polarized lasers and a magnetic field to lift 
the degeneracy in the sub-levels nip 42j. Our proposal simply extends this implementation 



through the addition of a third linearly polarized laser to facilitate |0) < — > \u). 

For performing a unitary transformation of BEC amplitudes in a pair of lattice sites, 
we utilize a scheme for the spin(state)-dependent transport of neutral atoms in an optical 



lattice 



27 



28] . By setting the wavelength Ai a ttice = 785nm, internal states |0) and |1) 
experience different corresponding dipole potentials Vo(x,9) = jV + (x,9) + ^V_(x,9) and 
Vi(x, 6) = V + (x, 9), where V±(x, 9) = V max cos 2 (kx ± 9/2), k = 27r/Ai at ti C o is the wave vector 
of the laser light propagating in the x direction, and 9 is the relative polarization angle 
between the pair of counter-propagating lasers. Hence for an atom in the superposition 
state a|0) + j3\l), increasing the polarization angle 9 will lead to a split in the spatial wave 
packet of the atom as it perceives a relative motion between the two potentials, resembling 
that of a pair of conveyor belts moving in opposite directions, each carrying one of the 
components a and (3. The relative displacement is given by Ax = #Ai att i ce /27r. 

Let us take a BEC initially prepared in the internal state |0) and distributed between 
two lattice sites \j,k) and \j,k') such that |^o) = &k\j, k) <g> |0) + ctk'\j, k') ® |0). We can 
now manipulate the amplitudes at and according to any desired unitary transformation 
in five steps depicted in Fig. [3h- (1) Using the three-photon STIRAP we apply a 7r-rotation 
to the BEC at \j, k) which transfers it entirely to the internal state |1) and the new state of 
the system becomes |"0i) = a.k\j, k) <S> |1) k') <g> |0). (2) Making use of the spin(state)- 

dependant transport, we increase the polarization angle by 9 = 2£(k — A;') 71 "/ Aiatticc causing 
the two wave packets to fully overlap at \ j, k') (selected as the stationary reference frame) and 
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hence \ip2) = \j, k') <8> («fc|l) + c*fc'|0)). (3) Using another three-photon STIRAP we perform 
an arbitrary unitary rotation R, this time at \j, k'), such that = \j, k') <S> (<5fc|l) + 5j^|0)). 
(4) Reversing the change in the polarization angle we transport the new BEC amplitudes 
&k and ay back to their original sites, i.e. |^4) = ctk\j, k) ® |1) + k') ® |0). (5) Finally 
performing another 7r-rotation on the state \j, k) we transfer the BEC back to the internal 
state |0) producing the desired outcome |^s) = &k\j,k) ® |0) + <Sfc'|j, k') <g> |0). Note that 
internal states |0) and |1) are only used to facilitate the pair-wise interactions and both BEC 
wave packets will be in their internal ground state |0) before and after they interact. 

This scheme can be readily extended to simultaneously activate all the pair-wise inter- 
actions required for performing the unitary rotations in Eq. [TJ We emphasize that all the 

(cj,) operations along the rows (columns) of the optical lattice are performed concur- 
rently, since the structure of the CS decomposition (Eq. [2]) is identical for all coin operators 
and changing the polarization angle 9 triggers the same spin (state)-dependent transport 
across the entire optical lattice. The effect of using different coin operators for each node 
appears in step (3), where the control STIRAP can perform different unitary rotations at 
various lattice sites. At the conclusion of the walk, BEC densities throughout the lattice can 
be detained via scanning electron necroscopy Q,Q or spin-selective absorption imag - 



mg 



431 ] . although the latter requires repeated runs of the experiment for each node density 



measurement. The corresponding quantum walk distribution is then derived by integrating 
the BEC amplitudes over an area A^tice x lattice centered around the key lattice sites. 
This will effectively include in the distribution, any residual amplitudes in the neighboring 
intermediate sites, which are nonetheless substantially lower than the amplitudes in key 
lattice sites and would therefore have a minimal effect on the final result. 

The proposed quantum walk scheme offers a polynomial speedup over an equivalent quan- 
tum circuit implementation, highlighting the expected trade off between resource and time 
scalability. A quantum circuit can in principle represent the walk's Hilbert space using 
m = log 2 (A/' 2 ) entangled qubits, which is by far more resource efficient. Then, implementing 
a generalized M 2 x M 2 unitary operator % Ci for each step of the quantum walk amounts to 



performing a m-qubit gate operation that can be realized with around 4 m CNOT gates [44]. 
Since the quantum circuit can perform at most m/2 simultaneous CNOT operations at any 
one time, each step of the quantum walk requires at least (4 m ) / ("72) = lM A j log2(A/" 2 ) oper- 
ational stages. This is compared to only M — 1 operational stages needed for implementing 
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Eq. [Q 

Spin(state)-dependant BEC systems have also been considered as serious contenders for 
building a quantum computer 45j . This is despite the acute sensitivity of the BEC internal 
states |0) and |1) to the external magnetic-field environment, leading to phase decoherence 
times that are presently in the order of a few ms 29(. Nonetheless, comparing this with 



a single-site transport time (~ 50/is) [271, |28| and STIRAP pulse durations (~ 60/is) [42], 



and also noting the successful realization of spin(state)-dependant BEC transport for up 
to 7 sites reported in j^j, a "proof of principle" implementation (i.e. the first few steps 
of the walk on an arbitrary graph with a few nodes) should indeed be possible, utilizing 
the existing experimental techniques. Since our proposed implementation scheme is in fact 
not inherently bound to any one physical system, naturally as this and other prospective 
quantum computing hardware grow in scale and fidelity of operations, so will the complexity 
of graphs on which the quantum walk can be performed. 



* Electronic address: wang@physics.uwa.edu.au 
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FIGURES 




FIG. 1: (a) A complete 6-graph. Any generalized graph can be constructed by removing edges 
(dotted lines) from the complete graph; (b) Quantum walk Hilbert space and a particular mapping 
^ni> k) — ► \k,j); (c) T is replaced by alternating the direction in which C is applied in successive 
steps of the walk. 




FIG. 2: Schematic diagram of a three-photon STIRAP operation in a 87 Rb atom. Internal levels |0), 
|1), \a) and \u) are coupled by three laser pulses with frequencies fii, and and polarizations 
that are linear, left circular cr_ and right circular <t+ respectively. 
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FIG. 3: (a) Steps for applying a unitary transformation to BEC amplitudes trapped in a pair of 
non-neighboring optical lattice sites; (b) The first three steps on a 2D optical lattice with BEC site 
separation 1 = 2. 
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